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We investigate theoretically the feasibility of an experimental test of the Bell-type inequality 
derived by Mermin for correlated spins larger than 1/2. Using the Schwinger representation, we link 
the output fields of two two-mode squeezers in order to create correlated effective spins between two 
observers. Spin measurements will be performed by photon-number-resolved photodetection, which 
has recently come of age. We examine the effect of nonideal detection quantum efficiency — and 
any other optical loss — on the violation margin of Mermin's inequality. We find that experimental 
violation is well accessible for spins larger than 1 for quantum efficiencies compatible with the current 
state of the art. 



INTRODUCTION 

The violation of Bell inequalities [ ] by measurements on quantum correlated systems with spins larger than the 
original 1 /2 value is of interest as a test of a new quantum regime, and possibly as a test of the quantum to classical 
transition [ ']. In 1980, David Mermin showed that the "standard" Bell inequality cannot be violated for spins larger 
than 1/2 and derived an alternative inequality which is violated by quantum mechanics for higher spins [.-,]. More 
recently, Collins, Gisin, Linden, Massar, and Popescu (CGLMP) discovered a generic method to investigate locality 
violations for higher-dimensional systems [.\\ — Qdits rather than Qbits [4] — that may also offer increased violation 
and resistance to loss. 

Classic spin-1/2 Bell inequalities have been violated under a variety of conditions over the past five decades. In 
particular, important experiments have closed various loopholes, such as the locality loophole [ - ], the fair sampling 
loophole [ ], and, most recently, the freedom-of-choice loophole [ ]. 

In this paper, our motivation is to pave the way for an experimental test of Bell-type inequalities for spins larger than 
1/2 in general. An experimental violation of Bell inequalities for spin-1 systems was demonstrated [1] but no other 
experimental work has been conducted on this subject to our knowledge. Here, we study an experimental protocol 
that was proposed earlier by Reid and Drummond [11, 12], albeit in the context of dichotomic (Qbit) measurements. A 
related but different proposal was also made by Gerry and Albert [ I '>] . This protocol (detailed in Section II) consists in 
using two pairs of entangled optical fields shared between two observers and viewed, in the Schwinger representation 
[14], as 2 spin s particles in the zero-spin state used by Mermin's high-dimensional Bell inequality. Here, we are 
interested in the experimental feasibility of violating Mermin's inequality for s > 1 and focus on a realistic model 
that includes arbitrary photodetection losses. This model is computationally intensive, hence we will focus solely on 
Mermin's inequality, which is the simplest computationally, even though an investigation of the CGLMP inequality 
(among others) is an equally interesting goal for future work. 

An example 

Despite the relative simplicity of Bell's original inequality, it resists straightforward generalization to higher spins, 
as was pointed out by Mermin [ ']. A good illustration of the situation is provided by examining the Clauser-Horne- 
Shimony-Holt (CHSH) inequality [ ], which can be constructed by taking the results of the measurements, by Alice 
and Bob respectively, of the components of two spins Sa and along axes denoted by angles a, j3, 7, and 5 about 
the z axis 

C = \mA,a + ■mA,^\mB,i3 + \mA,a - mA,f\ mB,5 (1) 

For spins 1/2, and normalizing the measurement results to values of ±1, it is easy to show that 

C < 2. (2) 

However, as is well known, this counterfactual CHSH inequality is violated by quantum mechanics, which allows 
Cmax = 2-v/2 for select angles when the spins are in a maximally entangled state. Let's now consider the general case 
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of 2 spins s. It is relatively easy to show that, in this case, a counterfactual inequality is 

C ^ 2s^ (3) 
Like in the spin 1/2 case, we now consider the quantum state of total spin 

l(ss)00) = -^i= ^ (-1)— |s,m)Js,-m),, (4) 

for which the spin correlation matrix elements take the simple form [2] 

{SA,aSB.ji) = -^s(s + l)cos(a-/3). (5) 
In this case the CHSH inequality becomes 

Icosfa — /3) + cos(7 — /3) + cos(a — 5) — cos(7 — (5)1 < . (6) 

s + 1 

The absolute value on the left-hand side of Eq. (6) is bounded from above by 2^/2, which corresponds to the maximum 
violation that quantum mechanics can offer. In that case, the inequality becomes 

O^^Om (7) 

which is clearly satisfied by s > 1/2. These spin values therefore pass the CHSH locality test. 

However, this simple example should only be regarded as the tree that masks the forest. Indeed, there are many 
Bell- type inequalities that are violated by quantum mechanics for higher spins. Reference [10] gives an instance for 
s = 1, along with an experimental violation, and the aforementioned Refs. [ , ] do so for arbitrary s. Moreover, 
work has also been done on Bell inequalities involving higher-dimensional systems in terms of dichotomic-variable 
measurements (e.g. photon-counting thresholds, photon-number parity) [11, 12, 17, 18]. 



Mermin's inequality [2] 

We now recall Mermin's counterfactual inequality, which is violated by quantum mechanics in the zero eigenstate 
of the total spin: 

)av (8) 

where 

s{\mA,a - mB,i3\)av = ^ |m - m' |P(to, to' , Q; - /3) . (9) 

m,m' 



In the quantum mechanical derivation, the zero-spin state of Eq. (4) yields 
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and, as aforementioned in Eq. (5), 

{mA^amB,'f)av = (5'a,q5'_b,^) = -7^s{s + 1) cos(q! - /3), (11) 

In evaluating the violation of Eq. (8), Mermin computed (numerical) optimal measurement angles for maximal viola- 
tion. The formulas above assume perfect measurements. In this work, we generalized Mermin's inequality to include 
the effect of measurement inefficiency in a particular optical implementation. The paper is organized as follows. 
In Section 2, we detail our experimental implementation proposal. In Section 3, we introduce our loss decoherence 
(detection quantum efficiency ry) model and apply it to Mermin's inequality. In Section 4, we present the results of 
our numerical evaluations for s = 1/2 to 9/2 and rj = 0.4 to 1. We then conclude. 
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PROPOSED EXPERIMENTAL IMPLEMENTATION 
Einstein-Podolsky-Rosen (EPR) entanglement between Alice and Bob 

We consider the experimental setup sketched in Fig. 1. Two two-mode squeezers, or optical parametric amplifiers 
0PA#1 and 0PA#2, are realized by using two optical parametric oscillators (OPOs) below threshold. Each OPA 
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FIG. 1. Proposed experimental setup for the violation of Mermin's inequality. OPA: optical parametric amplifier. PBS: 
polarizing beamsplitter. HWP: halfwave plate. PNRD: photon-number-resolving detector. Note that all optical paths in air 
between two consecutive PBS must be in phase, except for that of 02 which must be exactly out of phase. 

emits a two-mode squeezed, or EPR, state and each entangled mode is sent to each partner, Alice and Bob. Alice and 
Bob therefore each receive two independent fields, and each of their two fields is entangled with one of their partner's 
fields. This double EPR state is [19] 



n,=0 ^ 712=0 ^ 



(12) 



where ri_2 are the squeezing parameters. We now show how this state can be recast as a perfectly entangled spin 
state, even for finite values of ri 2- This was first pointed out by Drummond [ ] and Reid et al. [i-!]. 



Spin entanglement between Alice and Bob 



We use the Schwinger representation [ ' '], which consists in defining the following effective spin operators 

Sax = 2 (ai«2 + ai4) 
1 

Yi 
1 

2 



Sav 

Saz 



J. {a\a2 - ai4) 
1 ( t 



a\a2 



1 



(13) 
(14) 
(15) 

(16) 
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for Alice, and 



b2bl 



Sbx = ^ (4^1 

SBy = ^ (4^1 - b2b\) 

Sbz = I (4&2 - b\h 



Si 



b\bi+blb2 b%+hlb2 



(17) 
(18) 
(19) 

(20) 



for Bob. It is easy to check that these operators satisfy the canonical commutation relations of a quantum angular 
momentum. The usual eigenstates are those of (S*^, 5*2): 



for Alice, and 



I Sa,TOa), 



SB,mB )j 



nAi + nA2 nAi ~ nA2 
2 ' 2 



nBi + nB2 nBi - ^52 



(21) 



(22) 



for Bob. They can be measured in a straightforward manner if one can count photons with unity efficiency. The 
detectors in Fig. 1 will be photon-number-resolving superconducting transition-edge sensors, which have been installed 
in our laboratory at the University of Virginia, in collaboration with the groups of Sae Woo Nam at the National 
Institute of Standards and Technology and Aaron Miller at Albion College. Such detectors have been shown to allow 
photon-number resolution at quantum efficiencies exceeding 95% ['JU]. 

It is worthwhile to make explicit the physical meaning of these effective spin observables. Unsurprisingly, the total 
rotational energies given by observables g coincide with the total energy in Alice's and Bob's fields, respectively. 
The spin projection along z coincides with the photon number difference between the two fields on one partner's side. 
Finally, it is straightforward to see that the physical meaning of the spin projections along x and y is simply the phase 
difference between the two fields on one partner's side, since the observables S^^y are clearly interference terms. These 
considerations therefore completely specify the measurement procedure [21], which we detail in the next subsection. 

We can thus rewrite the double EPR state as 



s—O rn— — s 



(tanhri)"+'" (tanhra)''" 
cosh ri cosh r2 



\s m)j^\s -m)g. 



(23) 



This state can be put in an interesting form provided the squeezing parameters are identical ri — r2 — r. We then 
obtain 



EPR^ ) = ^ 



s=0 



(tanhr)' 
coshr 



2 s 



E 



IB ■ 



(24) 



This state is clearly a maximally entangled spin state at any given value of s, which can be heralded or defined 
by the measurement of the total photon number on Alice's and Bob's sides. In addition, we can obtain the exact 
zero-total-spin eigenstate by phase shifting either the A2 or the B2 optical path by tt: 



\n2)A,={-ir\n2) 



A2 ' 



which yields 



EPR^ ) = V2s + 1 



(tanh r)' 
coshr 



(ss)OO) 



(25) 



Thus, the double EPR state can be made equivalent to the eigenstate of two spins s coupled into a zero total spin. 
Do note that the maximal entanglement of the state is solely guaranteed by ri r2 = r, and in no way by the value 
of r, which does not have, in particular, to be large. 
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Spin measurement procedure 

In order to rotate the "spin analyzers" on Alice's and Bob's sides, we need to be able to measure the arbitrary spin 
component defined as 

Sa — COS a Sz + sin a Sx ■ (26) 

This is the purpose of the halfwave plates and subsequent polarizing beam splitters (PBS) in Fig. 1, which implement 
variable beamsplitters right before each set of detectors and can therefore rotate the measurement of the effective spin 
component from an intensity difference {S^) measurement to an interference (say, S^) one. Indeed, if both waveplates 
are at the same angle a/4 with respect to the PBS, the fields at each of Alice's detectors will be 

. a , a 

Ai = cos ^ «! + sm - 02 (27) 

. .a a 

A2 — ~ sm — oi + cos — a2 (28) 

and the expressions of the photon number on each of Alice's detectors will be 

a I ■ 2 t sin a / | ^ 



Nai = cos - ajai + sin — 0^02 H — [ala2 + aialj (29) 

. 2 Q! t 2 t sina / + +\ 
Na2 = sm — a[ai + cos — 0302 — [a[a2 + aiCjJ (30) 

so that the sum and difference measurements give 

Na+ = a\ai + 4a2 ^ Sa (31) 

Na- = 2 cos a Saz + 2 sin a Sax = '^Sa^o, (32) 

and the same of course on Bob's side. Hence, this experimental setup realizes the spin measurements of Eq. (26) for 
testing Mermin's inequality. 

We now turn to the core of this paper, i.e. the investigation of nonideal spin measurements, as defined by the 
nonideal photon-counting which stems from optical losses in the experiment. 

TESTING MERMIN'S INEQUALITY IN THE PRESENCE OF LOSSES 

Modeling loss decoherence 

We model loss decoherence by considering the effect of beam splitters of transmissivity r^j, j S [1,4], on the 
detected crosscorrelations. These beamsplitters can be assumed to be just before each detector but their position 
in the experiment between the source and the detector doesn't actually matter. They just serve as — valid — models 
[22, 23] for compounded propagation and detection losses. 

The situation is complicated by the fact that decoherence turns the pure state | <?!)) into a statistical mixture. We 
must therefore use a density operator description from now on. The statistical mixture is obtained from the density 
operator, partially traced over the loss ports of the beamsplitter: 

= Trjoss ports C^BS I 0) (01 C^BS ' (33) 

where Ubs is the unitary operator describing the action of 4 beam splitters, each with two input ports (one of which 
occupied by a vacuum field) and two output ports [one of which is traced over in Eq. (33)]. We now turn to the 
detailed expression of p' . 

Decohered Fock state 

Modeling the effect of losses as beamsplitters of transmissivity ry, we can use the well-known result [_ .] that a Fock 
state interfering with vacuum at a beamsplitter is turned into the "binomially entangled" output 

" / \ h 

UBs\n)jQ),,=Y.r]vHl-v)'^\k),\n-k),, (34) 



fc=0 



6 



In the case of decoherence, only one output port is available and the other must be traced over. The initial pure-state 
mode a is therefore turned into a statistical mixture described by 

P,-Trfc, [uBs{v)\n)jQ),,{QlA'^\aUlsiri)\ • (35) 
Evaluating this expression yields the familiar binomial probability distribution for the output photon number 



(36) 



Decohered arbitrary single-mode state 
We now consider any superposition of Fock states | ^ ) = V'n I ) • The initial density operator is thus 

\^P){^A = Y.'^"^n'\n){n'\ (37) 

n.n' 

and we can show, following the exact same procedure as before, that the partial trace associated with the decoherence 
process turns the operator \n) {n' \ into the operator 



" " / \ 5 / '\ 5 / / , 



E 



fc—max(0, n—n') 



k J \k + n' ~ n 



(38) 
(39) 



an expression which we can make more symmetric by using the sum index change k = k + [n! — n)/2 and the fact 
that max[±(rt — n')] = \n — n'\: 



lir. 



2 

E 



n 

11 — n' 



n 

n' — n 



(40) 



Decohered arbitrary two-mode state 
We now consider an arbitrary two-mode state 

IV) (^1= i'nun,rn'^,n'Jni){n[\(8>\n2){n'^\ 

ni ,712 -'"'1 ,^^2 



(41) 



The action of decoherence is independent on each mode, therefore we can directly use the result of Eq. (40) to show 
that the operator | ni ) ( n'^ | (g) | rt2 ) ( I becomes 



n. 



E 



Kl + -4 — / 



Ki + 



2 \ '"1 ' 2 



E 



K2 



"2 

712— 



H2 



%^(l-'?2)-^^+^ 



K2 + 



2 \ '^Z + ^ 



(42) 
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Now that we have two field modes, we can make use of the Schwinger representation to understand the efi'ects of 
optical losses on the effective spins that will be used in the final experiment. We find: 

E V'.,„^V'.V|sm)(s'm'|, (43) 

with the same definitions as Eq. (21) for the spin numbers. In addition, we replace the summation indices Kj in 
Eq. (42) with the spin variable (cr, n) such that 

- o + + ^2 + (44) 



M = n + o - ^2 - , ' (45) 



2 V 2 2 

1 
2 

These expressions define both kets in Eq. (42) as the spin ket \(t,^) and it is straightforward to show that the bras 
become the spin bra {a ^ s' — s, fi + m! — m\. Other helpful relations are 

s' — S + to' — TO , , 

= (T + M + ^ (46) 

s' — S — to' + TO , , 

K2=(T-^l+ . (47) 

The upper bound for cr is s and the lower bound depends on the relative numbers of photons ni and 712, as given in 
Table 1. 



TABLE I. Lower bounds cr,„i„ in Eq. (48). 

ni ^ n2 and n'l ^ n'2 s — s' 
ni ^ n2 and n'l ^ n'2 
m ^ na and n[ ni, 
m ^ nz and n[ ^ n'2 



Then, the action of decoherence turns operator | s m ) ( s' to' | in Eq. (43) into the operator 

1 11 ,1 

s + to\ ^/ s' + to' \^ f'^ ^ ^\ ^ ( s' — m! 

(T + fj, J \<j + fi + s' — s + m' — m J \a — IJ, J \a ~ fi + s' ~ s ~ ni' + m 
I cr, /i ) ( CT + s' — s, /i + m' — TO I (48) 



The expression in Eq. (48) is thus Alice's decohered | s to) (s' m' | operator. Note that all s,to, s',to' must be the 
same for Alice and Bob because of the quantum EPR correlations before decoherence (and, again, the finiteness of 
squeezing plays no role in this). This yields the final decohered entangled density operator, which will be used to 
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compute all averages in Mermin's inequality: 

CO CO S s' S (7a S (7f^ 

^'-EE«vE E Mr-""'-'"' E E E E 

S — Os'—O ra — — Sm' — — s' (7a—(T,nin^a — — (7a(7b—(7minfJ'b — — (7i, 

'lal 'Ia2 'lal) ['^ 'Ia2) 

i 11 

s + m \ ^ / s' + m' \ ^ / s — TO \ ^ / s' ~ m' 

O-a+f^aJ \(Ja + f^a + S' - S + rn' - mj \<Ta-fiaJ \Cra ~ f^a + s' - S - m' + 171 
I CTq, )^ ( 0-a + s' - S.^la+m' -m\ 

s + to\^/ s' + m' \^/s — to\^/ s' — m! 

ffe + Mfc/ + /ifc + S' - S + to' — TO/ \0-;, - ^b/ - + ^ 5 - m' + TO 

I (Tb, ^6 )5 ( (76 + s' - S, /ib + to' - TO I , (49) 

where r^^ = (tanhr)"/ coshr. 



Computing Mermin's inequality in the presence of losses 

In order to evaluate the inequality in Eq. (8), we derive the expressions of the dual-measurement probability 
(left-hand side term) and of the crosscorrelations (right-hand side term). 



Dual-measurement probability 

We now derive the probability distribution of measurement results {sa, ma) and (sb, TOb), with respective "analyzer" 
directions a and /3, in order to evaluate the left-hand term of Eq. (8). Note that each measured "spin" Sa.b may 
independently be an integer or half an odd integer. The probability distribution is 

P{Sa,ma,Sb,mi,;a,/3) 

= Tr [p' {e-'^'^y \sama}{sama\ e^"^") ® (e"^'^^- | Sb ^6 ) ( Sb mb | e'''^-)] (50) 

= E-'.--'vE(-ir™"'-"' E E 

S,S' m,m' (Ta f^a <7b Mb 

'7al ^a2 ' Vbl Vb2 

(1 - 77al)""""^"+'^+"(l - 7?a2)"""+^"+'""(l " r]bl)""'"'"'+'+'" (1 - ,yb2)""''+^''+''"" 

1 11 .1 

s + TO \ ^ / s' + m' \ ^ / s — TO \ ^ / s' — m! 

<ya+ P-a) + ^ia + s' - S + - m/ VcTq - Ma/ V^a " Ma + s' - S - m' + TO 

s + to\^/ s'-|-to' \^/s — to\^/ s' — vn! 

c^b + fJ-bJ \(^b + Hb + s' — s + m' — mj \ab ~ ^J-bJ V^b — Mfc + — * — to' + TOy 

( CTa + s' - S, Ma + to' - m I e"'"'^» I Sa, Wq ) ( Sa, | e"^'^" I CTa, Ma ) 

{ab + s' - s, ^b + m' - m \ e"*^'^" | Sb, TOb ) ( Sb, TOb | e'-'^^y \ ab, ^J,b ) (51) 
We now use the angular momentum selection rules for rotation matrix elements, which yields 

aa + s' - S = (Ja = Sa (52) 

ab + s' - s ^ ab = Sb (53) 
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If s = s', this gives aa.b = Sa,b, respectively. However, if s 7^ s' , no suitable aa,b exists, which considerably simplifies 
the expression by introducing a Sg^s' term. Note that m' — m is always an integer. We obtain 

P{sa,ma,Sb,mb;a,f3) = 

E 4 / -|\2s-m-m' \ ^ Sa+A'o+ '^ 2~" ""2"*"" "* 2~" ""2"'""' 

~sr 2-^ \ ^) 2^ 'lal 'Ia2 'Ibl 'Ib2 



1 1 1 

s + m \ ^ I s + m' \^ I s — I s — m' 



Sb+tib+s — m 



1 11 

s + / s + m' \^/s — m\^/ s ~ m' 

Sb + fJ'bJ \Sb + fib + m' - mj \Sb-fibJ \Sb - fib - m' + m 

( Sa, Ha+m' -m\ e"*"^« \Sa,ma) { Sa, nia I e*"^" \ Sa,[la) 

( Sfe, ^6 + m' - m I e"*''^" | Sb, mt ) ( Sb, \ e'^^" | s,,, /if, ) . (54) 

This expression may already be evaluated, but for simplicity we will focus on the particular case where all losses are 
equal rjj ^ rj, Vj. In that case we get the slightly simpler expression 

P{sa,ma,Sb,mb;a,l3) = 

s m,m' fJ'a,f^b 

1 

s + m \ ^ / s + m' \ " / s ^ m \ ' / s — m 



Sa + Ma/ \Sa + Ma + m' - mj \Sa - fi-a J \Sa - fla - m' 
1 11 

s + m \ ^ / s + m' \^/s — m\^/ s ~ m' 



ySb + fJ-bJ \Sb + Hb + 1^' — mj \Sb — HbJ \Sb - fJ^b — m' + m 
with the well-known expression for the rotation matrix elements 



<:+m'-m.m„(«) («) d^+rW -rn.mAl^) dZ,mM, (55) 



4n,m2(") ^{smi\e "'■^y\sm2] 



/(s + mi)!(s-mi)! / / a\™i-™2 „ _^ ^ 

^V (. + m2)!(.-m2)! r^2) ^^2) ^''^ 

where the last symbol denotes a Jacobi polynomial. Note that rj = 1 forces 2s — Sa — Sb = and, in fact, s = Sa = Sb 
since the photon numbers and therefore, spin magnitudes, are perfectly correlated between Alice and Bob. This is 
of course not the case when 77 < 1 as photons are randomly lost between the various detectors. The ideal efhciency 
should of course give the same result as Mermin's original inequality, which wc will verify numerically in the next 
section. 



Cross correlations 



We now derive the expressions for the crosscorrelations on the right-hand side of Eq. (8). The arbitrary spin 
component of Eq. (26) can be written in terms of spin ladder operators 

S + S 

Sa = COS a Sz + sin a (57) 

which we can use in the crosscorrelation expression 



Cap = {SA,a Sb,I3) 



(58) 
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We then evaluate these matrix elements using the decohered entangled density operator p' . 

C^p = Tr [p' {e'°'^^ySA:,e~"^^^y) ® (e^'^'^^^S's^e-'^^^'')] 
= Tr 



p' [ cos a Saz 
cos/3 Sbz 



2 

sin (3 



{Sa+ + Sa-. 
{Sb+ + Sb- 



Again, we take all losses equal to simplify the computations. This yields 

OO S \ S ^ cL s 

c.p=Y.<r E (-if^^-™'cos«cos/3 E E E 



E 



s—Q ni- 



2((T„+CT6)('-| _ \2(-(T„-(T6+2s + 2m) 



s + m\ / s ~ m 



fa + PaJ \Cra — l^-a ) \Oh + Mb/ \Ob — Mb 



Ma Mb 



„ s a a s Ub 

I ^ -|^^2(s-m)-i Sinasmp E E E E 
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^a—^min ^^a—~^a ^b—^min ^b — ~^b 
2((T„+CT6 + 1)/-| _ \2(-(Ta-(T6 + 2s + 2m) 



s + m\ ^/ s + m + l\ ^/ s — m \^ j s — m — 1 

+ Ma/ \0'a + Ma + 1/ VCTa - Ma/ \^a - fJ'a - 

1 11 

s + to\^/s + to+1V^/s — m\^ / s — m— 1 
o-b + Mb/ + Mb + 1/ \o-b - Mb/ \o-b - Mb - 1 



v/caCcTa + 1) - Ma(MQ + 1) Vcrb(o-6 + 1) - Pb{^J■b + 1) 

^ -|^^2(s-m)+i Siiiasm/3 E E E E 

^2((T„+(T6-l)j--|^ _ ^>|2(-cr„-cr6+2s+2m) 

i i 1 

s + m\ ^/ s + m— l\^/s — 771 \ ^/s — m-f 

O-a + Ma/ \0-a + Pa - ^ J \(Ja ~ Ma/ \^a - Ma 
S + TTi \ ^ /■ S + TO — 1 \ ' I s — m 



<Jb + IJ-b 



<yb + ^-^b-^J \a-b - ^J■b J \a-b - Mb + 1 



s — m 



V^^aifTa + 1) - Ma(MQ - 'i-)V'^b{(Tb + 1) - fJ'bifJ'b - 1)| 



(59) 



(60) 



(61) 



We are now ready to make use of Eqs. (55,61) to evaluate Mermin's inequality quantum mechanically as a function 
of the loss parameter, or detection efficiency, rj. 



CALCULATION RESULTS 



We evaluated Eqs. (55,61) numerically for a range of detector angles (a, (3), efficiencies rj, and squeezing parameters 
r. After defining an angle a — 7r/2 = /3 — 7r/2 = 0as Mermin does, our expressions reproduce Mermin's results 
exactly at perfect efficiency r/ = 1. Evaluation of these probabilities at perfect efficiency is straightforward because 
the sum over variable s reduces to a single term, as the (1 — 7y)2(2s-s„-Sb) ^gj-j^^ gQgg ^ero for all but the single spin 
in question. In this regime, the experimental results are exactly equivalent to those achieved by Stern-Gerlach style 
measurements on a spin s particle. 

In the case r; < 1, the sum in the probability over variable s cannot be neglected. To evaluate these expressions, 
we cut off the sum over s at a point determined dynamically in the computation, depending on the squeezing 
and the efficiency. Then, we carried the sum out further to verify that our results converged to good precision. 
The numerical evaluation of these expressions placed heavy demands on computational time, especially for high- 
squeezing, low-efficiency, or high-spin systems, all of which required summing over many values for s. For these 
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difficult computations, we made use of the "Dogwood" parallel computing cluster managed by the University of 
Virginia Alliance for Computational Science and Engineering. Most of these expressions were evaluated on 1-4 
computing nodes, each equipped with two four-core, 3 GHz Intel Xeon EMT processors with 3 GB of RAM. Typical 
run time was anywhere from several minutes for the least demanding jobs to over 150 hours for the most demanding. 

Influence of decoherence on spin analyzer angles 

Figures 2 and 3 show the violation of Mermin's inequality, which is defined as the difference of the right-hand 
and left-hand sides of Eq. (8), for different efficiencies 77 at fixed spin s and OPO coefficient r. The inequality is 
violated when the "Violation" is positive in Figs. 2 & 3. The efficiencies decrease from the top to the bottom of the 
plots — higher efficiencies violate over a greater range of angles, and higher spins require higher efficiencies to produce 
violations. 



Violation 




FIG. 2. Plot of violation versus angle for different efficiencies for the s = 1 case. Higher values constitute greater violation; 
negative results can be explained by a local hidden variable model. Decreasing efficiencies produce decreasing violation. Here, 
r = 0.5. 



Violation 




FIG. 3. Plot of violation versus angle for different efficiencies for the s = 2 (four photon) case. The results are much the same 
as in Fig. 5, except the magnitude of violation is lower. Note that there is also a smaller angular range of violation and a 
greater sensitivity to inefficiencies. Again, r — 0.5. 

Unsurprisingly, a lower efficiency 77 results in considerably reduced violation margins for Mermin's inequality; as 
detailed below, violation also decreases with an increase in either the squeezing parameter r or the spin s. Additionally, 
the optimal angle for violation appears to shift very slightly as the efficiency decreases. 

It's also important to recall that the situation considered here is perfectly symmetric. All of the detectors have 
equal efficiency, and the choice of optimal angle 6 (defined by Mermin [2] ) depends crucially on the equivalence of the 
two detectors: the otherwise independent detector angles a and /3 are identical in this model. Further optimizing the 
violation of Mermin's inequality over the (really four-dimensional) parameter space of the angles might yield better 
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results for a particular choice of experimental conditions. In the results that follow, however, we only evaluated the 
inequality at the optimal angle for perfect detection. Thus, the violations we predicted here are, at worst, lower 
bounds for actual violations. 



Influence of decoherence on violation for different spins s 

After choosing to evaluate these results at the optimal angles for rj — I, we can examine the effect of inefficiencies 
more directly by plotting the violation against the efficiency, for r — 0.4, Fig. 4, and r = 0.2, Fig. 5. These plots 
are constructed by taking a single point (the maximum in the case of perfect efficiency, near the maximum for lower 
efficiencies) and combining these points in a plot of violation vs. rj. Different lines on the same plot represent, from 
top to bottom, different spins from s = 1/2 to s = 9/2 in steps of 1/2. 



Violation 




FIG. 4. Plot of violation versus efficiency for different spins s. Again, higher values constitute greater violation. Decreasing 
efficiencies produce decreasing violation, and lower spins violate more strongly than higher spins. Here, r = 0.4. 



Violation 
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FIG. 5. The same plot as Fig. 4 for r = 0.2. Note that the general trends are the same, but violation is stronger. Generally, 
higher values for r give lower violation for r\ <\ and vice versa. 

As Figs. 4 & 5 show, as the efficiency decreases, all but the lowest spins fail to violate the inequality; the higher 
the spins, the more deleterious the effect of inefficiency on violation. In Fig. 5, though, all spins tend to violate more 
strongly than in Fig. 4, especially at low efficiencies. This trend is general — the lower the squeezing parameter r, the 
less effect inefficiency seems to have on violating the inequality. 
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Joint influences of decoherence and squeezing 

Many of our results can be summarized in the plot of Fig. 6, which again shows the dependence of the violation 
of Mermin's inequality on r and rj. Note that Figs. 4 & 5 are cross-sections of the Fig. 6 at constant r, and that 
Figs. 2 & 3 have their maxima represented as single points on Fig. 6. Crucially, violations of Mermin's inequality are 
predicted for all of these spins under a range of experimentally attainable conditions. 




FIG. 6. This plot compares the effect of both squeezing r and efficiency i] on the extent of violation. The nine surfaces show 
the predicted violation for the nine different spins. The top surface is that for s = |- The next surface is that for s = 1, and 
so on, through s = 9/2. Note that the surfaces are labeled with the same color scheme as in Fig. 5 and Fig. 4 — these two plots 
are just the cross sections through the two planes r — 0.2 and r = 0.4 in this figure. As a reference, a plane is also plotted at 
zero violation. 

There are several trends in these data. First, the extent of violation of the inequality decreases monotonically as 
the detection efficiency decreases. This was expected. The extent of the decrease, however, depends both on the 
amount of squeezing and on the total spin itself. On the one hand, low spins almost always violate under the range 
of conditions plotted here, and have little dependence on either parameter {r] or r). Higher spin systems violate 
less strongly — even in the ideal case of perfect efficiency — and the extent of violation decreases much more rapidly 
with r] than for lower spin systems. That exacerbated sensitivity of high-spin quantum systems to decoherence is not 
too surprising either. However, it is important to mention here that several inequalities since Mermin's time make 
violation more robust by remedying one or more of these issues. Garg and Mermin later produced an inequality whose 
angular range of violation does not decrease with increasing spin [25], and more recent work has produced simple 
inequalities that are violated even for arbitrarily large spin [3, 26, 27] approaching the maximum possible violation for 
a 2 X 2 X d Bell- type experiment [28-30]. The CGLMP inequality [3] in particular, produces large violations even for 
large spins and also appears to be resistant to noise in the form of Werner states. We've begun exploring the CGLMP 
inequality in our optical framework, but the form of the inequality makes incorporating loss more computationally 
intensive than in Mermin's 1980 inequality. In future work, we plan to address the optical realization of these newer 
high-spin inequalities. 

Finally, we see from Fig. 6 that a given spin s produced by a strongly squeezed source is much more sensitive 
to decoherence than the same spin s produced by a weakly squeezed source. In the case of perfect detectors, the 
amount of squeezing doesn't affect the extent of the violations, but a weakly squeezed state actually produces greater 
violations than a highly squeezed state if the efficiency is low. This is a much less intuitive result if we associate a more 
highly squeezed state with less classical behavior, and thus expect a more highly squeezed state to produce greater 
violations. Our result can be simply explained by the fact that the photon-number distribution of the output fields 
depends strongly on the squeezing parameter: in EPR states as produced by the OPO, the probability amplitude of a 
given Fock state is Tg^r = (tanhr)"/ coshr , which in both limits goes as a decaying exponential in r. So, if the OPO 
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is barely squeezing, the probability for generating a total photon number, or spin state s, is much larger than the 
probability of generating even the next larger spin, whereas, as r increases, all of the photon-number states become 
equally probable, as displayed in Fig. 7. Inefficiencies thus effect a sort of "crossover" between the photon-number 



FIG. 7. Plot of the relative probability of producing a certain photon number in a Fock state as a function of the squeezing r. 
As the squeezing increases, generating larger and larger photon number states becomes more and more probable. 

states, with a strength increasing as the inefficiency increases: the beamsplitters that model losses annihilate photons 
from the signal field, which reduces the total spin corresponding to that Fock state in the Schwinger representation. 
Then, because the measured results are partitioned into individual spin-s subspaces before checking for violation, the 
measured distribution of projections in a given subspace depends both on the distribution in that subspace and on the 
distribution in higher subspaces, because the higher subspaces will couple to the lower subspaces with an interaction 
strength determined by the beamsplitter reflection coefficient, or inefficiency of the detector. The relative amplitudes 
of the two subspaces are therefore crucially important: on the one hand, even if there is strong coupling between 
the two subspaces through an inefficient detector, if the lower spin subspace has a significantly higher associated 
probability amplitude, then the effect of the higher spin subspace will be negligible, and this is what happens in the 
low-squeezing case. If, on the other hand, there is high squeezing, then the two spin subspaces have comparable 
probability amplitudes and there is a large effect if the detectors are inefficient. For imperfect detectors, higher 
squeezing thus causes many of the same effects as lower efficiency because both cause a lower-spin subspace to have 
increased coupling to a higher spin subspace. 

Nonetheless, it's important to remember that there is a significant practical difference between increased squeezing 
and decreased efficiency. Although theoretically, lower squeezing produces higher violations given inefficient detectors, 
with low squeezing the OPA produces high photon number states much less frequently. So, there's a tradeoff: low 
squeezing will produce higher violations than high squeezing in a nonideal detection case, but one will have to wait 
much longer to see the violations in an experiment. Increasing the efficiency allows for increased squeezing, which is 
crucially important for violating any high-spin inequality in a reasonable experimental timeframe. 



Our results show that it appears quite feasible to violate Mermin's high-spin Bell inequality under a range of 
experimentally viable conditions. After converting the photon number distribution produced by an optical parametric 
oscillator from the Fock basis to the spin basis, we find that the state could be written as combination of singlet states 
typically used to violate Bell-type inequalities. Incorporating inefficiencies into our probabilities by considering the 
effect of loss through beamsplitters allows us to determine the effect of actual experimental parameters on the outcome 
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CONCLUSION 
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of our tests. In particular, we find that Mermin's inequality is violated fors = l/2tos = 9/2 under a wide range of 
experimentally realizable efficiencies and squeezing parameters. 

In the future, we hope to build on these results by performing an experimental violation of these inequalities using 
the setup described previously. We hope also to examine other high-dimensional inequalities such as that of Collins, 
Gisin, Linden, Massar, and Popescu[ '] which — though more computationally intensive in this framework — may offer 
greater violations and increased resistance to experimental imperfections like loss and noise. 

This work was supported by NSF grants PHY-0960047 and PHY-0855632, and by the FINESSE program of 
Raytheon BBN, under the DARPA "Information in a Photon" program. RE was partially supported by an Echols 
Scholar Ingrassia Research grant at the University of Virginia. 
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